How to do THD analysis in MATLAB